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Abstract 

We derive the homogenization limit for time harmonic Maxwell’s equa¬ 
tions in a periodic geometry with periodicity length rj > 0. The considered 
meta-material has a singular sub-structure: the permittivity coefficient in 
the inclusions scales like and a part of the substructure (corresponding 
to wires in the related experiments) occupies only a volume fraction of order 
the fact that the wires are connected across the periodicity cells leads 
to contributions in the effective system. In the limit ry —>■ 0, we obtain a 
standard Maxwell system with a frequency dependent effective permeability 
and a frequency independent effective permittivity Our formu¬ 
las for these coefficients show that both coefficients can have a negative real 
part, the meta-material can act like a negative index material. The magnetic 
activity yr®® / 1 is obtained through dielectric resonances as in previous pub¬ 
lications. The wires are thin enough to be magnetically invisible, but, due 
to their connectedness property, they contribute to the effective permittivity. 

This contribution can be negative due to a negative permittivity in the wires. 

Keywords: Maxwell’s equations, negative index material, homogenization 
MSC: 78M40, 35B27, 35B34 

1 Introduction 

Light and other electromagnetic waves exhibit interesting refraction phenomena in 
meta-materials. While no natural material posesses a negative index, by now, there 
exists a variety of negative index meta-materials (a meta-material is an assembly 
of ordinary materials, arranged in small substructures). A smart choice of the 
microscopic geometry can lead to macroscopic properties of the meta-material that 
are not shared by any of the ordinary materials that it is made of. Even though 
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no natural negative index material is known, a distribution of metallic or dielectric 
resonators and metallic wires can act, effectively, as a medium with a negative index. 
We present here a periodic meta-material (periodicity length t] > 0) with singular 
sub-structures (wires with relative radius of order r], i.e. with absolute radius of 
order 77 ^) which is, on the one hand, close to the experimental set-up, and, on the 
other hand, accessible to a rigorous mathematical analysis. 

From a modelling point of view, the Maxwell system is very simple: We in¬ 
vestigate solutions {E'^, H^) to the time harmonic Maxwell’s equations in three 
dimensions. The system contains three parameters: the frequency a; G M and two 
material parameters, the permeability fi = /i(a;) and the permittivity e = e:(x). 
Since natural materials have a relative permeability close to one, we assume that 
the permeability coincides with that of vacuum, /r = /iq > 0. All the complex 
behavior of the micro-structure is encoded in one single coefficient, the relative per¬ 
mittivity — )■ C. Here, we consider periodic coefficients Sr^ with large values 

of \er^\ in the inclusions. Using the parameter £0 > 0 (permittivity of vacuum), we 
investigate solutions (U^, of 

curlU^ = , (1.1) 

curlif^ = —ioje^EQE'^. (1.2) 

Our aim is to analyze the behavior of the solutions {E'^, H'^) in the limit 77 —)■ 0. 

The relative permittivity is given by two complex numbers and the geometry 
of two types of inclusions, S,, C and T^ C M^, both periodic with periodicity 
77 . The set Hrj consists of bulk inclusions, the number of which is of order 77 “^ (one 
inclusion in each periodicity cell). The set T^ represents a system of long and thin 
wires, each wire has a radius of order 77 ^ and length of order 1 (the wire has a 
macroscopic length and a width that is small compared to the periodicity length). 
With parameters eb,Sw G C with A(efe), A(e^) > 0 we study the coefficient 

{ £fe 77“2 in , 

£^ 77 “^ inT^, (1.3) 

1 in \ (S^ U T^) . 

Our result is a description of the weak limit {E, H) of a solution sequence 
(U^, H'^). We show that an appropriately dehned pair (U, H) := (U, [jl)~^H) solves 
a Maxwell system with two effective parameters /i®®(a;) and cp. fll.bp - fll.7p . 
The two effective parameters are given by explicit formulas that involve cell prob¬ 
lems for the electric and for the magnetic held. In an appropriate geometry and for 
an appropriate frequency w, we hnd that both coefficients can have a negative real 
part (simultaneously): The effective permeability can become negative due 

to resonances in the bulk inclusions Due to their vanishing volume fraction, 
the wires do not inhuence the parameter 77 *^®(a;). On the other hand, due to their 
special topology (they form connected objects of macroscopic length), they inhu¬ 
ence the permittivity If Ew has a negative real part (which is the case for many 
metals), then can be negative, see fll.lip . 
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1.1 Literature 

Half a century ago, Veselago investigated in the theoretical study ca materials with 
negative e and negative ja. Since the product Sja is positive, waves can travel in such 
a medium, but surprising effects such as negative refraction and perfect lensing can 
be expected. Since no natural materials exhibit a negative /i, the studies of Veselage 
have not been continued until in about 2000 hrst ideas were published on how to 
construct a negative index meta-material, see e.g. [16]. Regarding applications of 
negative index materials we mention the effect of cloaking by anomalous localized 
resonance, see e.g. [B] and [H]. 

A mathematical analysis of meta-materials became possible with the develop¬ 
ment of the method of homogenization. The homogenization technique was suc¬ 
cessfully applied in many situations, ranging from porous media to wave equations. 
Also the homogenization of Maxwell’s equations has been performed. The standard 
result of a homogenization process is the following: Given periodically oscillating 
coefficients ar^{x) and a corresponding solution sequence m^(x), every weak limit 
u of the solution sequence satishes the original equation with an averaged coef- 
hcient a®®. In the context of the Maxwell system, the oscillating coefficients are 
arf{x) = {en{x), firtix)) and the solution is In [T8|, results of this 

kind are obtained for Maxwell’s equations. 

Of particular interest are those homogenization results that lead to a qualita¬ 
tively different equation for the limit u. Examples are the double porosity model in 
porous media [2] or the dispersive limit equation for waves in heterogeneous media 
[9l [Tn|. The problem at hand is similar in that we want to combine positive in¬ 
dex materials to obtain an effective negative index material. Typically, such effects 
are obtained by micro-structures that involve extreme parameter values and/or by 
micro-structures that contain finer substructures. 

In this work, we will actually use both, extreme values and fine substructures, to 
obtain a surprising new limit formula in the homogenization of Maxwell’s equations. 
One of the pioneers in the held is Bouchitte who, together with co-authors, initiated 
the held with the analysis of wire structures 0. Ha- It was shown that extreme 
coefficient values in the (thin) wires can lead to the ehect of a negative 

In order to obtain a negative index material, additionally a negative p®® has to be 
created. Several ideas have been analyzed. Based on the analysis of a model prob¬ 
lem, a mathematical result has been obtained in [13]. A truely three-dimensional 
analysis of a setting that is also close the some of the original designs has been 
carried out in [7]. In that work, it was shown that the periodic split ring structure 
can lead to a negative /i®®. A similar result has been obtained later also for hat 
rings of arbitrary shape in [TB] (while the rings in [7] had to be tori). 

While this seems to be the state of the art with regard to metallic inclusions, a 
simpler approach using dielectric materials has been invented in [3] and later devel¬ 
oped in [3]: Dielectric inclusions can lead to resonances (so-called Mie-resonances) 
which result in an ehective magnetic activity. With resonant dielectric inclusions, 
a negative p®® can be obtained even without subscale variations of the periodic 
geometry. 
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With this contribution we close a gap that has been left open by the above 
mentioned works: The emphasis has always been to create either a negative 
mm or a negative p®® [31111 El [15] — here we present a construction that achieves, 
simultaneously, a negative e®® and a negative /r®® (we always refer to the real part 
of the coefficient). An important point in our construction is the decoupling: e®® 
depends only on and the shape Sy of the bulk inclusions; given Sy, we can choose 

to have £®®^ negative (independent of the frequency). Tuning the frequency, we 
can generate a resonance in Sy which makes /i®® negative; this process does not 
affect £®®. 

We conclude this section with a comparison of our results to those of [S], where 
also the effect of thin wires has been analyzed. In [5], the volume fraction of 
the wires is much smaller than in our study, namely 6^^ ~ exp(—2p“^) (here, it is 
9r^ ~ we observe that the permittivity is scaled as in our work as \er,\ ~ 9^^)- The 
effective equation in [5] is a Maxwell type system that is coupled to a macroscopic 
equation for a new variable (denoted by J in [5]). Our result is much simpler in the 
sense that we obtain a standard Maxwell system in the limit, see fll.bp - fll.7p . In 
this sense, our result is very different to that of [5], despite the similarities in the 
setting (thin wires with vanishing volume fraction of macroscopic length). In order 
to understand how such a different outcome is possible, let us mention the following 
point (somewhat technical, but important): the tiny radii in [5] make it impossible 
to construct test-functions as in fl3.15p - fl3.17p . hence the macroscopic limit cannot 
be obtained as in our setting. 

1.2 Geometry 

The underlying idea of our approach is simple: We use dielectrical bulk-inclusions 
(given by the subset S^) as Mie-resonators, they lead to a negative /i®®. Addition¬ 
ally, we include thin wires (given by the subset T^) that have a negative £ (which 
is the case for many metals); the wires lead to a negative £®®. Since the wires are 
thin, their effect decouples from the Mie-resonance. 

We emphasize that, dehning the two sets and T^, the Maxwell system is 
completely described by fll.ip - fll.3l) . Below, we consider the limit of any bounded 
solution sequence and do not specify boundary conditions. Our result is therefore 
applicable to any boundary value problem. 

Microscopic geometry. We construct a periodic meta-material in three space 
dimensions. We start from the periodicity cube Y = [0,1)^; since we will always 
impose periodicity conditions on the cube Y, we may also regard it as the flat 
torus T^. We now construct two subsets, Sy,ry C Y, compare Figured] for an 
illustration. 

Dielectric resonator. The dielectric resonator is given by a subset Sy C Y. We 
assume that Sy is a simply connected open set that has a Lipschitz boundary and 
that is compactly contained in (0,1)^. We call the set the dielectric resonator since, 
in order to have resonances in Sy, we must assume that Sb has a positive real part, 
CP. flTTOD. 
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Thin wires. The second subset are the (metallic) wires. We assume that three 
wires T’lf GY,j = 1,2,3, are contained in Y. We assume that each wire has a 
radius of order r] inside Y and that it connects two opposite sides of dY in a periodic 
fashion. We assume additionally that the wires do not intersect each other or the 
resonator Sy; this assumption is not satished in some of the experimental designs 
(“hshnet-structure”), but we do not see it as an essential property for our method 
to work. 

In order to make our assumptions precise (and, at the same time, to keep the 
construction of special test-functions in Section |3] accessible), we restrict ourself to 
cylindrical wires of relative radius a > 0: For a point G (0,1)^, the wire 

with index j = 3 has the central line Fy := x [0,1) and is given by 

Ff :=i?„,((|/f\#))x[0,l)cF. (1.4) 

The wires with indices 1 and 2 are constructed accordingly. The union of the three 
wires is denoted by Ty '■= Uj=iry'^- The construction assures that the relative 
radius of the wires inside the cube is ay such that the wires vanish in the limit 
?7 —)■ 0. We will see that the wires do not enter the cell problem, but, due to their 
connectedness across cells, they do affect the macroscopic equations. The real part 
of the permittivity in the wires has to be negative in order to obtain a negative 
quantity in fll.lip . Since metals can have a negative permittivity, we think of 
metallic wires. 




Figure 1: Left: The macroscopic domain and the scatterer R G Ll (here, 
is a ball and i? is a prism). Maxwell’s equations are solved on fl. The equations 
contain the coefficient which is 1 in \ i? and which has large values and a 
micro-structure in R. Our interest is to describe the solution helds {E’^, H'^) in the 
limit 7] — )■ 0. Right: The periodicity cell Y with the subsets Sy (the resonator) and 
Fy (the wires). In the two subsets the coefficient has absolute values of order 
7 ]“^. The wires have the radius ay in the periodicity cell, which means that they 
have the radius ay^ in the macroscopic domain. The periodic construction yields 
connected wires with a length of order y^ = 1. 
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Macroscopic geometry. We study Maxwell’s equations in an open set C 
Contained in is a second domain R G Q with R G Q. The set R consists of 
the meta-material, on the set fl \ R we have relative permeability and relative 
permittivity equal to unity, see Figure [H left part. 

In order to dehne the microstructure in R, we use indices k E 1? and shifted 
small cubes Yj] := T]{k + Y). We denote by /C := {/c G G R} the set of 

indices k such that the small cube Y^ is contained in R. Here and in the following, 
in summations or unions over k, the index k takes all values in the index set /C. 
The number of relevant indices has the order |/C| = 0{t]~^). 

Using the local subsets Sy and Ty-^, we can now dehne the meta-material by 
setting 

S, := U + S'-) ■ r, := U U + r«) =\Jq(k + rp. ( 1 , 5 ) 

keic keic je{i,2,3} keic 


1.3 Main result 

With the help of cell-problems, we dehne two tensors G The tensor 

e®® is dehned in fll.lip and has two contributions: The matrix A®® is dehned with 
the standard electrical cell problem, the additive term aem can make the real part 
of e:®® negative. The ehective permeability tensor /i®® is dehned in fl2.22p through 
a coupled cell-problem that appeared already in EE] and was studied there. The 
formula for p®®(a;) is a result of the resonance properties of Sy; a negative real part 
of /i®®(a;) can be expected for frequencies that are close to eigenvalues. 

Effective equations. We derive an effective Maxwell system in H: 

curl E = iufiojlH, (1.6) 

curl H = —iue^sE . (1.7) 

In this system, the coefficients are fi{x) = e{x) = 1 for x & Q\R, whereas for x E R 

i{x) = , jl{x) = . 

A consequence of our result is the following: Since outside the scatterer R G fl the 
ehective coefficients are equal to 1, the held {E,H) coincides there with the weak 
limit {E, H). In particular, for x E kl\R, as r] ^ 0, the solutions {E^, 77^) converge 
to {E,H), solution to the ehective system fll.6^ - 01.7p . Let us now formulate the 
precise statement in our main theorem. 

Theorem 1.1. Let {E^,H^) be a sequence of solutions to fll.ip - fll.2p . where the 
coefficient Sn is given by (na and the geometry of inclusions S,, and wires T^ is as 
described above, cf. dn and mai, We assume that the solution sequence satisfies 
the energy-hound 

[ \Hfi^ + \e,\\Efi^<C, 

Jn 


( 1 . 8 ) 
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where C does not depend on rj, and that converges weakly in L‘^{Q) to a 

limit {E, H). Then the limit {E, H) has the property that the modified fields E := E 
and H := fl~^H solve fll.6p - fll.7p in the sense of distributions on 

1. Comments on the fields E and H. The field E coincides with E, we intro¬ 
duced the new name only to have a consistent notation. Loosely speaking, the field 
H is the limit of the fields outside the inclusions. It is obtained from the two- 
scale limit Ho{x,y) of by geometric averaging, see 02.25^ : geometric averaging 
is introduced in Section 12.31 

2. Easy parts of the Theorem. Since we assume that the imaginary parts of Sb 

and are positive, the estimate fll.Sp implies the boundedness of E’^ and 

Since every weak limit is also a distributional limit, we may therefore take the 
distributional limit in (HU and obtain curl E = iup^oH. Because of E := E and 
H = jlH, we conclude fll.bp in the sense of distributions. The challenging part of 
the proof is to derive the effective equation fll.7p from the Maxwell equation fll.21) . 

3. Comments on the a priori estimate fll.Sp . We have formulated our main 
theorem with a boundedness assumption on the solution sequence. This is done for 
convenience, the discussion of the a priori estimate is available in the literature. We 
give here only a brief summary of some facts. 

(a) It is sufficient to make only the (weaker) assumption of L^(r2)-boundedness 
of the solution sequence, 

[ + \Efi^ < C. (1.9) 

Jn 

Indeed, assuming only fll.91) and TsSb, TsEu, > 0, one obtains fll.Sp on subdomains 
n with R G Cl CC by testing the equations with the solutions (multiplied with 
cut-off functions) and by considering real and imaginary parts of the results. 

(b) Analysis of the scattering problem. In a scattering problem, one considers 
n = and uses an outgoing wave condition for the diffracted field as a boundary 
condition. In this setting, it is not easy to derive a local L^-bound as in fll.9p . In [7] 
and [15] we succeeded to work with a contradiction argument: Assuming that the 
local L^-norms t^^ are unbounded along a solution sequence, the rescaled solution 
sequence {E’^, H’^) /t^j is locally bounded; the homogenization result for bounded 
sequences can be applied and provides a limit system. In the case f,, —)■ oo, the 
limit system has a vanishing incident field and therefore only the trivial solution. 
Together with a compactness argument, this provides a contradiction to the fact 
that the solutions {E’^,H‘’^)/tr^ are normalized. 

(c) We mention that the trick described in (b) did not work in the case of obsta¬ 
cles as in [8], where the inclusions have a macroscopic length. The long obstacles 
in that contribution result in the geometrical problem that two points in R with 
distance of order rj can have a distance of order 1 in the metric that is induced by 
the coefficient 8^- In the contribution at hand, the complement of the inclusions is 
connected in each periodicity cell, which suggests that the scattering problem could 
be solved with Theorem O arguing as in [7]. 

Interpretation of the resnlt and comparison to experiments. The result 
of this contribution is the following: One can combine in each periodicity cell a 
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resonator and a metallic wire strnctnre. The resonator with relative permittivity Sb 
can lead to a negative magnetic coefficient A negative relative permittivity 

Eiu in the wires can lead to a negative electric coefficient Since the wires are thin, 
the effects are decoupled: The coefficient is determined by cell-problems that 

do not see the wires, the coefficient is a result of a negative permittivity in the 
wires and is independent of the frequency. 

In experiments, the metallic wires can indeed have a large negative real part, 
hence our assumptions on the wires are well justified. On the other hand, the experi¬ 
ments typically use metallic resonators (e.g. split ring resonators); in this set-up, the 
resonance is not based on the Mie-resonance of dielectrics, but on capacitor-inductor 
interactions. The methods in this article can also cover this case: Combining the 
thin wires with the metallic split-ring resonators of [7] or [TJ] would lead to the 
same result (we emphasize that the coefficient k in [7] or [12] can also be a complex 
number). We present here the case of dielectrics, since the analysis of the resonance 
is much simpler in this case. 

Let us high-light why both coefficients can have eigenvalues with negative real 
part. The tensor /r®® of fl2.22p can be expressed with eigenfunctions of an 

eigenvalue problem in the cell as 

^ S ' -fy ("‘' X 

where k = and (A„)„eN are the eigenvalues. The formula is taken from [3], 

it is identical to the present case, since the cell problem is identical. We sketch the 
essential arguments that lead from fl2.22p to fll.lOl) in Appendix The formula 
shows that, if Sbk'^ is close to an eigenvalue of the cell problem, then the coefficient 
/i®® can be large in absolute value and the real part can have both signs. 

For the effective permittivity £®®, we have the formula 

:= A®®-1-, (1-11) 

where the (positive) tensor A®® is defined in 03.231) . While A®® (and, obviously, also 
vra^) are always positive, a negative sign of Ree^ can lead to a negative real part 
of (all the eigenvalues of) e:®®. 

Notation. Spaces of periodic functions are denoted with the sharp-symbol jj, 
e.g. as H^(Y). We use the wedge-symbol for the cross-product, v A w := {V 2 W 3 — 
V 3 W 2 , V 3 W 1 — V 1 W 3 , V 1 W 2 — V 2 W 1 ), and the curl-operator curln = V A n. Constants C 
are always independent of 77 , they may change from line to line. 

2 Two-scale limits and cell problems 

In this section, we define the (rescaled) dielectric current and consider the two- 
scale limits EQ{x,y), HQ{x,y), and Jo{x,y) of the three sequences E^, and J^. 
We derive cell problems that determine the two-scale limits up to the macroscopic 
averages E{x) and H{x). The procedure uses two-scale convergence and follows 
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closely the lines of other contributions in the held, e.g. PEIE]. The crucial dif¬ 
ference regards the inclusion of thin wires. In our setting, the cell problems are 
identical to those of [3] — despite the wire structures. On the one hand, this fact 
allows to use the results of |3] on the cell problems. On the other hand, the effect 
of the wires decouples from the effect of the dielectric inclusions. The wires enter 
only in the derivation of the macroscopic equation. 

The flux J^. Besides and we consider a third quantity, namely the 
rescaled dielectric held : O —)■ C^, deflned by 

J’^:=r)e^E\ ( 2 . 1 ) 

The deflnition of and estimate fll.Sp imply 

( 2 . 2 ) 

Jn JnxcEnur^) JEr, V Jrr, V 

The L^(0)-boundedness of the unknowns and implies that we can And 

a sequence rj = rji ^ 0 and two-scale limit functions Eq{x, y), Hq{x, y), and Jo{x, y) 
such that E’^, and converge weakly in two scales to the corresponding limit 
functions (compactness with respect to weak two-scale convergence). 

2.1 Cell-problem for Eq 

Lemma 2.1 (Cell-problem for Eq). Let {E’^,EI^) he a sequence of solutions of the 
Maxwell’s equations as described in Theorem 1 1.11 and let EQ{x,y) be a two-seale 
limit of E’^. Then, for almost every x E R, the funetion Eq = Eq{x,.) G H^(X) 
satisfies 


curly Eq = 0 in Y , (2.3) 

divy Eq = 0 in Y \ Sy , (2.4) 

Eq = 0 in Ey . (2.5) 

For a given cell-average E{x) = fyEQ(x,y) dy, there exists a unique solution Eq of 
the above problem. 

Easy parts of the proof and consequences. We claim that there exist three ele¬ 
mentary solutions E^,E‘^,E^ G H^{Y) with the canonical normalization, 

EE,Y —)■ , E^ solves fl2.3^ - 02.5p and J E^ = ej . (2.6) 

The elementary solution E^ can be constructed with a solution of a Laplace problem: 

Let 0-^ be harmonic in T\Sy with Q^{y) = -y-ej in Ey. Setting E^{y) := VyQ^-\-ej 
provides a solution with average = Cj. Since every solution is given by a 

gradient of a potential, the construction implies also the uniqueness of solutions to 
prescribed averages. 
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As a consequence of the lemma, we obtain that the two-scale limit EQ{x,y) can 
be written with the help of the three basis functions as 

3 

Eo{x,y) = J2Ej{^)E^iy)^ (2.7) 

i=i 

where Ej coincides with the j-th component of the weak limit of the sequence E'^. 


Proof of Lemma \2.1\. It remains to derive the cell problem and to show that it is 
not affected by the presence of the wire structure. The properties of two-scale 
convergence imply rjcvxlE'^ —)• cvxly EQ{x,y) while, on the other hand. Maxwell’s 
equation fll.ll) and boundedness of yield r/curl 77’' —)■ 0. This provides fl2.3p . 

When restricted to the function ri~^E^ is L^-bounded by estimate fll.Sp . 
This implies fl2.5p . 

Outside U T^, the electric held E^ has a vanishing divergence. This implies 

divj^ i7o = 0 in Y \ (Sy U Ty) (2.8) 

for the one-dimensional wire-centers Ty = ljj=i ^y in the sense of distributions. 
We claim that fl2.8p implies fl2.4p . In order to prove this claim we have to argue 
with the fact that the set Ty has a vanishing i7^(y)-capacity. 

We work here only with p = 2, i.e. with the 2-capacity, for dehnitions see e.g. m- 
The one-dimensional hbres in dimension n = 3 have a vanishing capacity. This can 
be concluded from Theorem 3 in m, page 254, since n — p = 3 — 2 = 1 and 
the one-dimensional Hausdorff-measure of the wires is hnite. The more elementary 
argument is that points in two space dimensions have a vanishing capacity and the 
corresponding cut-off functions can be extended to three-dimensions for straight 

By the dehnition of a vanishing capacity, for every h > 0, there exists a function 
Ts with the properties 

: r ^ R is of class H^{Y) , (2.9) 

vanishes in a neighborhood of Ty , (2-10) 

||1 “ T('<5||hi(v) < ^ • (2T1) 

Let now (p G C^{Y) be an arbitrary test-function with support contained in 
Y \ Ey. We consider (^T^, for which fl2.8p and the properties of imply 

0 = [ Eo{y)-V{pTs)= [ Eo{y)-VpTs + 0{6)^ [ Eo{y)-Vp (2.12) 


>Y 


>Y 


'Y 


as 5 —)■ 0. We hence obtain fl2.4p . 


□ 


2.2 Cell-problem for (i^o, Jo) 

The cell-problem for the magnetic held is obtained as in the case without wires, 
we do not have to use vanishing capacity arcuments. In fact, the cell problems are 
identical to those in the case without wireqj. 

'^Our equation (12.151) coincides with equation (14) of [3], except for the negative sign in front of 
the second term; we correct here a typo in (14) of [3]. We remark that the analysis of [5] proceeds 
with the correct sign after equation (14). 
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Lemma 2.2 (Cell-problem for the pair Let be a sequence of 

solutions as in Theorem M . 1\ and let be defined by fl2.ip . Let HQ{x,y) and Jo{x,y) 
be two-scale limits of and Then Hq{x, .) satisfies, for a.e. x E R, 


curly 77o 4- iujsqJo = 0 m F, 

(2.13) 

divy 77o = 0 m F, 

(2.14) 

while Jo{x ,.) satisfies, for a.e. x E R, 


curly Jq — iSbOJyioHQ = 0 in Sy, 

(2.15) 

divy Jo = 0 m F, 

(2.16) 

Jo = 0 inY \ Sy. 

(2.17) 


Proof. The Maxwell equations fll.ll) - fll.2p immediately imply relations fl2.13p - 
fl2.16lh one only has to use the properties of two-scale convergence and the dehnition 
of J’’^. We note that these hve equations hold either on all of Y or on Sy, hence 
the justihcation is as in the problem without wires. 

Relation fl2.17l) is a consequence of the a priori estimate: On the set : = 

Q \ U r^), the function ri~^J’^ = = E’^ is bounded. More precisely: 

is bounded in L‘^{Q) and hence converges strongly to 0 in L^(0). 

This implies for the two-scale limit function the relation 

Jo{x, y) = 0 for almost every y eY \ (Sy U Ty) , (2.18) 

for almost every x E R. This relation is identical to fl2.17p . since Ty is a set of 
vanishing measure. □ 

It was shown in [3] that, with an appropriate normalization, the cell-problem 
fi 2 .i 3 p - fi 2 .i 7 p has a unique solution. A natural normalization is to prescribe the geo¬ 
metric average ^ 77 of the cell solution, we dehne this average in the next subsection. 
The normalized solutions are denoted by = H^{y), j = 1, 2, 3, 

Hfi F ^ C^ {Hfi R) solves flOTH - flOTll and ^ W = Cj . (2.19) 

As a result, the two-scale limit Hq{x, .) of can be written as a linear combination 
of the three shape functions Hfi Denoting the coefficients by Hj{x), we write 

3 

Ho{x,y) = Y,Hfix)W{y). ( 2 . 20 ) 

It remains to establish the relation between the weak limit H (x) and the coefficients 
77, (x). As the weak limit of 77^, the function 77(x) coincides with the F-average of 
77o(x,.). Taking the F-average of 02.2011 . we obtain 

77(x)= [ Hfix,y)dy = y^Hfix) [ W (y) dy = H (x) (2.21) 

Jy .1 Jy 
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by setting 

j HHy)dy. (2.22) 

Relation fl2.2ip was used in Theorem II.II to define the field H. 


2.3 Geometric averaging 

There are two possibilities to average a function u : Y ^ C^. The standard 
averaging procedure is the volumetric average, given by the integral over Y. In 
a geometric averaging procedure, we seek for the typical value of an integral over 
a curve that connects opposite sides. The concept has been made precise with 
the definition of a circulation vector in [3]: To a vector field m G X = {m G 
H^(Y ; C^)|curlyM = 0 on F \ Sy}, we associate the circulation vector / m G C^. It 
is defined by the property 



(2.23) 


for every test-function ip G L^(F, C^) with (p{y) = 0 for almost every y G Sy and 
with V • = 0 in the sense of distributions on Y, more precisely: Jy p ■ Vf/ = 0 for 

every U G H^{Y). 

We must verify that fl2.23p indeed defines a vector f u G C^. To this end we 
must show that, for every function p with average JyP = Ci, the expression JyU-p 
has the same value. By linearity, it is sufficient to show that fyU-po vanishes for 
every test-function pq with a vanishing average. 

As a preparation, we note that the conditions on the test-functions p imply that, 
for every function [/ G II^(Y), there holds fy^^y^ p ■ VU = ^yP ■ Vt/ = 0. Since 
Y \ Sy is simply connected, a curl-free function u can be written, up to an additive 
constant uq G C^, as a gradient of a periodic function U, i.e. u{y) = uq + VU{y) 
ioT y eY \ Ey. We therefore find fyU-po = /y\y('Wo + VU) ■ po = uq ■ fy po = 0. 

To understand the definition of f u G through fl2.23p better, let m be a 
function in X and let 7 be a curve that does not touch Ey and that connects the 
lower and the upper face of Y, i.e. the side {y^, = 0} with {y^ = 1}, in points that 
are identified by periodicity. If u is continuous and 7 is differentiable, there holds, 
for the 3-rd component. 



where r is the tangential vector along 7 . The formula expresses that every compo¬ 
nent of the circulation vector of u coincides with the corresponding line integral of 
u. Due to curlyM = 0 in F \ Ey, the line integral is independent of the curve 7 by 
the Stokes’ theorem (we recall that F \ Ey is simply connected). 

Let us demonstrate fl2.24p with the help of special test-functions. The test 
functions are constructed with the help of cylinders that are similar to the wires 
in our homogenization setting. Focussing on the third component, we use T^ := 
Bs{{yi,y2)) X [0,1). Choosing the point {yi,y2) G (0,1)^ sufficiently close to the 
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boundary and 5 > 0 small, we accliieve that does not touch Sy. In order 
to evaluate the circulation vector, we consider the function (p(y) = eslpgiy) that 
points in the third coordinate direction and vanishes outside the cylinder. The 
function (p has a vanishing divergence and it vanishes on Sy, it is therefore a valid 
test-function. Its volume average is Jyip = Vol 3 (r 5 )e 3 =: VsCs. We can therefore 
calculate the left hand side of fl2.24l) as 


63 • 




-/ 

^ J Bs({yi,y2)) 


u 


63 dz^ 


'^(21.22) 


dzi dz2 


U ■ T . 


In the hrst equality, we used the dehnition fl2.23p . in the second equality, we used 
Fubini’s theorem and decomposed the integral over the cylinder into an integral over 
the cross-section and integrals over the corresponding lines. In the third equality, 
we used that all line integrals coincide and that V 3 = Vol 3 (r 5 ) = 1- We 

have thus obtained 02.241) . 


Application to the two-scale limits Eo and Ho. We can now understand bet¬ 
ter the connection 02.211) between H{x) and H{x). Confronting the normalization 
of H^{y) in 02.19p with the dehnition of /i®® in 02.22p . 


/ 




= e 


3 ! 



(2.25) 


we see that the matrix /i®® corresponds to the factor between Waverages and geo¬ 
metric averages of cell solutions. Instead, for the cell-solutions ^ the two averaging 
procedures coincide. 


/ 


E^ 



= e 


3 5 


(2.26) 


which is a consequence of curl^i?-^ = 0 on F: all line integrals coincide, even if the 
curves intersect Sy. This fact yields E{x) = E{x). 


3 Macroscopic equations 

In this section, we obtain the macroscopic equation 01.71) . The starting point for 
the derivation is the Maxwell equation 01.21) and our knowledge about the two-scale 
limits Eq and Hq. The main idea is to construct special test-functions for 01.2p . A 
hrst type of test-functions has a vanishing curl in the whole domain, a second type 
of test-functions has a specihed curl in one of the wires. 

3.1 The ? 7 -dependent test-functions avoiding wire integrals 

We start the construction by dehning, for j E {1, 2, 3} hxed, a potential : V —)■ M 
of class 0^ E H^^(V). We set := F\(I]y Ty*) and dehne 0^ as the solution 
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of 


0 ^( 2 /) =-ej-2/in Sy uljr^*, (3.1) 

0 ^( 2 /) = 0 in , (3.2) 

AQi^iy) = 0 in Y; . (3.3) 

The construction allows to define —)■ R”' with the gradient of 0^, 

< = vej + ej . (3.4) 

This provides a function that satishes 

curl'd^ = 0 in y, (3.5) 

divd:^ = 0 in Y *, (3.6) 

d'^ = OinEyUlJr^', (3.7) 

i¥=j 

= ej in Ty^ , (3.8) 


with the normalization Jy = e^. We note that is closely related to the solution 
■. Y ^ R” of the electrical cell problem (which does not see the wires). In fact, 
we have the following convergence result. 

Lemma 3.1. There holds 


d)) —>■ in L‘^(Y) as rj ^ 0. (3.9) 

Proof. We fix the index j G {1, 2, 3} and use variational arguments. The potentials 
0^ minimize the Dirichlet functional 

A(0) := J |V0p on Xn := {0,, G Trj|(y, R)|0,, satisfies 03.1 p and 03.21) } . 

The potential 0-^ (that provides the cell solution E^ = + ef) minimizes the 

Dirichlet functional A on the set X* := |0 G i7}(y, R)|0 = —ej ■ y in Sy}. From 
Xy C X* we immediately obtain A(0-^) < A{Q{^). 

Our aim is to show the convergence of the energies as —)■ 0, 

A(0);) ^ A(0^'). (3.10) 

We observe that 03.101) implies the lemma: The functions 0^ are bounded in 
H^{Y), hence we can select a weakly convergent subsequence. The weak lower 
semi-continuity of A on X* together with 03.101) implies that every weak limit of 0)) 
minimizes A, it therefore coincides with the (unique) minimizer 0T Since O3.10p im¬ 
plies additionally the convergence of norms, we can conclude the strong convergence 
of the sequence. By uniqueness of the limit, the whole sequence converges. 
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The convergence fl3.10p follows if we show, for arbitrary > 0, that the relation 
A(0^) < A(Q^) + ^0 holds for every sufficiently small rj. In order to obtain this 
inequality, we construct, for arbitrary 5o > 0, a comparison function 0^ G with 
- 4 ( 8 ') < - 4 ( 8 ') < ^( 8 ') + S„. 

We use the cut-off functions of fl2.9p - fl2.lip . more precisely: The function 
is periodic, vanishes on the subset of all wires Ty = ULi equals 1 on Sy and 
is close to 1 in the iffj|(y)-norm. We additionally use a function T| with similar 
properties: is periodic, vanishes on the subset Ty* of two wires, is equal to 

1 on Ty-^ U Sy and close to 1 in the iifj|(y)-norm. We construct 

eJ-e'T, + (-e,-!/)(l-TJ). (3.11) 

The functions 0;^ are periodic on F; we emphasize that this fact is only true since 
we do not set 0^ = —yj on the j-th wire. Furthermore, the comparison functions 
satisfy 0^ G for rj > 0 sufficiently small. They have the Dirichlet energy 

.4(8',) := ^ I ve'/ = ^ I V(e' T,) - V(s,(l - Ta)|2. 

The smallness assumption fl2.1ip and the boundedness of 0-^ (due to the maximum 
principle) imply 


Aiei) < |V(0^' T,)|2 + 0{5) < |V0^f |T,|2 + 0{5) < A(0^) + 0{5 ). 

Because of the minimality of 0)) we have A(0))) < A(0^) and obtain therefore 
A(0^) < A(0'^) -|-0((5). This yields fl3.10l) and thus the statement of the lemma. □ 

Proposition 3.2. Let be a sequence of solutions as in Theorem \l.l[ Let 

(f G be arbitrary and let be the special test-functions defined above. 

Then there holds, as rj ^ 0, 


'n 


H''^{x) ■ {dUx/rj) AVifi^x)) dx ^ / / HQ{x,y) ■ {E^[y) AVg){x)) dy dx , 


n Jy 


(3.12) 




enE’^{x)f}Ux/r])ip{x) dx ^ / / Eo{x,y) ■ E^(y) ip{x) dy dx . (3.13) 


n JY 


Proof. Relation fl3.12p follows from the strong two-scale convergence of ^{^{x/ri) to 
E^{y) that was shown in Lemma [3.11 For the convergence of products of strongly 
and weakly two-scale convergent sequences we refer to [T]. 

Concerning fl3.13ll we note that = I holds in \ (S^ U F,^) and that 

vanishes on (S^ U F^) \ F)). Again, the strong two-scale convergence of 'di^^x/rf) to 
E^{y) implies the convergence to the double integral. □ 
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3.2 The calculation of wire integrals 

With the last proposition, we are almost in the position to derive the macroscopic 


limit equations: The left hand side of fl3.12j) coincides with tested against 


a suitable test-function. The left hand side of fl3.13p coincides with tested 

against the same test-function. Since the two quantities are related by the Maxwell 
equation (O, a comparison of the two right hand sides can provide the missing 
effective equation that relates curl H with E. 

There is one missing piece: the left hand side of fl3.13p does not include the 
integration of Sr^E'^ ■ Cj over the j-th wire, T^. Its limit cannot be calculated as in 
the previous subsection, since T;^ connects opposite sides of the cell in a periodic 
way. Technically: = —ej -y cannot be demanded in fl3.1l) . since this function 

is not periodic in yj. 

We therefore calculate the missing integral with a different approach. The result 
is, in some sense, not very surprising: The wire integrals of E'^ converge to line 
integrals of the two-scale limit function Eq (up to the factor 7ra^ that measures the 
volume of the cross section of the wire). This stability of line integrals over E'^ in a 
two-scale limit process is a consequence of the fact that the curl of E^ is controlled. 

Proposition 3.3. Let {E'^,H'^) be a sequence of solutions as in Theorem \l.l\ and 
let ip G Cc (^j arbitrary. Then there holds, for j = 1, 2, 3, 



(3.14) 


Even though the result is suggestive, the proof of Proposition 13.31 requires some 
preparation. As in the proof of Proposition 13.21 we must construct special test- 
functions. 

Geometry and notation. We will construct an ? 7 -dependent sequence of test- 
functions Y —)■ in several steps. To simplify notation we assume the follow¬ 
ing: By choice of other coordinates, the unit cell is Y = (—|, |)^, the coordinates 
are gathered as y = (|/i,|/ 2 , 1 / 3 ) = iy^y^)- We consider j = 3, i.e. that wire inside 
Y that runs in ea-direction. We exploit that the wire is straight and assume for 
convenience that it has the central line ^ = 0, i.e. Py = {?/ G = 0} and 
Py^ = Bar^iVy) = i?or)(0) X [—\, |). We choose a radius 5 > 0 sufficiently small 
such that Bs{Ty) does not intersect the other sub-structures, Py^, Py^, and Sy. In 
the following, we only consider r] with ay < 6/2. 

Construction of test-functions. We start with a scalar function 


fjri ■ [ay, 00 ) -)■ M, := aylog{r)ps{r ), 


(3.15) 


where ps : [0, 00 ) —)■ M is a smooth and monotone cut-off function with ps{r) = 1 
for every r < 6/2 and p&{r) = 0 for every r > 6. The function allows to define 
the vector potential 0,,: 


9 , R, ^rfy) ■■= iPMul) <^3 ■ 


(3.16) 
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With the help of the potential 0^, we can hnally dehne the test-function by 
setting 

—curly©^ for \y\ > at], 


g^i :Y ^ 


9viy) ■= ' 


1 

ar] 


yi 

\ 0 / 


for \y\ < ag. 


(3.17) 


For another explanation of the ideas, let us describe the construction as follows: 
We consider a (truncated) two-dimensional fundamental solution outside the disc, 
and dehne gn outside the wire as its rotated gradient. Inside the wire, we dehne grj 
as a rigid rotation. 


Properties of the test-functions. We start with the observation that gn is 
continuous, coinciding with a normalized tangential vector on the cylinder sur¬ 
face. Indeed, using the two-dimensional tangential vector r = (— 1 / 2 ,2/i)/|^| = 
(-|/2,2 /i,0 )/|^|, we have dehned gn{y) = (- 2 / 2 , |/i, 0)/(ar7) = r for r = |^| = ag. On 
the other hand, for |^| > ag, we hnd with r = |^| 

9r)i.y) = -curlj^©^ = dr'ipni'r ')^, and dril^nir) = ag^psir) + ag\og{r)drP5 {r). 

In particular, because of drPs = 0 in a neighborhood of r = ag, there holds 9 n{y) —t 
T as r = \y\ ^ ag. 

The function Pn is bounded in the unit cube, Hfi'jjHoo < 2 independently of 0 < 
g < go for some small r^o > 0. This property follows by inspection of the expression 
for dr'ipni'^)^ which the hrst term is bounded by 1, the second term is bounded 
by Cg\ log(77)|. Furthermore, there holds 

||5'r?||L2(y) ^ 0 . (3.18) 

Indeed, for the integral over the small region |y| < ag, we can exploit boundedness 
of Pn- For the integral over \y\ > ag, we calculate for the hrst term j^{g/rY r dr = 
g‘^\og{5/g) —)■ 0 as ?7 —)■ 0. The second term is uniformly bounded by Cg\ log(? 7 )| 
and we conclude 03.1811 . 

In the next step, we calculate the curl of Pn- For \y\ < ag, there holds 

curly pniy) = — 63 = — eg . (3.19) 

ag V / “h 

Instead, for ag < |^| < 5/2, we hnd, using div©^ = 0, 

curln Pniy) = ^(^,^ 2 ) Hlog(l^l)] eg = 0. 

Finally, for |^| > 5/2, there holds 

^curly Pniy) = [«log(l^l)P^(l^l)] eg =: Gs(y). (3.20) 
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We extend Gs by setting Gs{y) = 0 for |^| < 6/2. We then find ri~^curly gr^{y) = 
Gs{y) for 1^1 > arj. It is the main point of our construction, that this function is 
non-singular in the limit rj ^ 0 , despite the factor ri~^. 

As a last step in this preparation, we observe that Gs{y) points always in direc¬ 
tion 63 and that it is independent of y^. Its average can be calculated with Stokes’ 
theorem, using the normal vector u = {yi,y 2 , 0 )/{ar]) on the cylinder surface and 
the fact that n A r = 63 : 


Gs{y) -e^dy = 


>Y 


/ -curly Qyiy) ■ 63 dy = 

' \y\>ar) V 


1 


/ ^driiy) ■ esdSiy) 

'\y\=ar] V 

(3.21) 


=ari V 


- dS{y) = —27ra . 


Convergence of wire averages. With the help of the above oscillatory test- 
functions we can now prove the convergence result. 


Proof of Proposition\3^ From fl3.19p . we know the curl of inside the wire. 
Transforming into the x-variables, we have curl^, [grtix/rj)] = 2 /[arf)e 3 in F^. We 
can therefore express the wire integral with the help of the curl of an approri- 
ately designed test-function and proceed with a straight-forward calculation. In 
the equation marked with “ fl3.20p ” below, we perform an integration by parts of 
the curl-operator, in the limit marked with “ fl3.18p ” we exploit the boundedness of 
curia;which follows from the hrst Maxwell equation and the L^-boundedness on 

— [ -^E'^{x) ■ 63 ip{x) dx [ E^{x) ■ curh [gyix/r])] ip{x) dx 
“ Jr^ V ir3 



ils:^ 


E^{x) ■ curia; [gy^x/r])] p{x) dx — 


curl^E'^{x) ■ [gy{x/ri)] ^p{x) dx 



E'^{x) ■ curia; [gy^x/v/)] (p{x) dx 



E'^{x) ■ [gy{x/ri)] A Va^p{x) dx 



arm. 


arm. 


' E^[x) ■ Gs{x/g) (p{x) dx 
Eo{x,y) ■ Gsiy) (p{x) dydx 


Eo{x,y)-e 3 ] j Gs{y) dy p{x) dx 


Q JY 


= - E 3 {x) / Gsiy) ■ e 3 dy(p{x)dx 


lY 


27ia [ E 3 {x) ip{x) dx . 


This calculation shows the limit fl3.14p and hence Proposition 13.31 


□ 
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3.3 Derivation of the macroscopic equations 

Limit process in fll.lj) . We can take the distributional limit of fll.ljl and obtain, 
in the limit rj ^ 0, 

curl E <— curl E’^ = —)■ iufioH = iu/aojlH. (3.22) 

We recall that the last equation is a consequence of the definition of fi and H. 
The above distributional limit equation already provides fll.61) . the first of the two 
effective equations. 


Re-writing of the two-scale limit integrals. It remains to conclude the sec¬ 
ond effective equation, cm. We will obtain this equation from fll.2p . exploiting 
Propositions 13.21 and 13.31 In order to prepare the calculation, we re-write terms 
that have been obtained in Proposition 13.21 
We define the coefficient matrix 


A^-.= j E‘{y)-E={y)dy. (3.23) 

With this definition, we can write the Wintegral on the right hand side of fl3.13p . 
for X E R, as 


^ E„{x, y) ■ E’{y) dy = ^ U EEiy) ■ E^iy) dy'j E.{x) = (A’^Eix)), . (3.24) 

To calculate the right hand side of 03.121) . we use the expansion 02.20p of Hq 
and the definition of the circulation vector: The function E^ A is a test-function 
which vanishes on Sy and which has a vanishing divergence; regarding the latter 
we recall E^ = V0-^ +ej which implies, for fc = 3, V ■ {E^{y) Ae^) = diE{ — d 2 E{ = 
did 2 Q^ — 92^1©'^ = 0. This allows to express the integral of a product with the 
circulation, 

/ Ho{x,y) ■ {E^{y) Aek)dy ^^ / H\y) ■ {E^{y) A eu) dy 

Jy Jy 

(3.25) 

Hi{x) ei ■ {cj A Ck ). 

With this preparation, we can now perform the limit process. 


t A(x) EHy) A = t 


Limit process in 01.2p . In order to perform the limit r; —>■ 0 in 01.21) . we use an 
oscillating test-function. We choose a smooth function <^9 : hi ^ M with compact 
support and fix j G {1, 2, 3}. We consider = 'diiix/rj) (p{x) withfrom 03.4p . 
Then the second Maxwell equation 01.21) yields 


curl ■ ipr) = —ioJSo / e^E’^ ■ 


(3.26) 
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It remains to evaluate the limits of both sides of fl3.26l) . We start with the left 
hand side. In the subsequent calculation we use hrst integration by parts, then 
curly= 0. In the limit process we exploit fl3.12p of Proposition 13.21 


curl ■ ip„ = / H'^ ■ curl = 


H^{x) ■ {'dUx/'t]) A V(p(a;)) dx 


m 





iih:^ 


HQ{x,y) ■ {E\y) ^V<^{x))dydx = 


H{x) ■ (cj A V(p(x)) dx 


n JY 


= / H{x) ■ curl {(p{x) Cj) dx= (curl.^) ■ ej cp . 


in 


We now calculate the right hand side of 03.261) . In the hrst equality, we use that 
'd{^{x/ri) vanishes on and on all P^ with i ^ j, and that it coincides with ej in 
P^. The limit process y ^ 0 for the two integrals has been prepared in 03.13P and 

mM- 


EriE’’^ • ijj^ 


ErjE^ ■'d{^{x/ 7 ])ip{x) dx + / ErjE'^{x) ■ ej(p{x) dx 


i ri 


iis:^ 



Jn\{T.r,urr,) 

Eo{x,y) ■ E^(y) p>{x) dydx + ira'^Eii; / Ej{x)p>{x)dx 


n JY 


in 


A®"E(x) • Ej <p{x) dx + Tra^Euj / E{x) ■ Cj <p{x) dx . 


Jn Jn 

Since j G {1, 2, 3} and tp = ip{x) are arbitrary, we obtain from 03.26p 

curl.^ = —iojEQ^A''^ + Tia^Eyj) E{x). 

This is exactly the effective equation 01.7p . since we have dehned the effective per¬ 
mittivity in 01.lip as e®® = A®® -|- tiq^Ew 


Conclusions 

We have investigated Maxwell’s equations in a periodic material with small period¬ 
icity length 7] > 0. The permeability is set to 1, the permittivity is assumed to have 
extreme values of order 0{r]~‘^) in the periodic inclusions, it is 1 outside the inclu¬ 
sions. Two types of inclusions are present: bulk inclusions and wire inclusions. The 
dielectric bulk inclusions have an impact on the effective permeability /i®®, an effect 
that has been studied before in [3]. In our setting, the cell-problems for /r®® are 
identical to those of [3] and the study of the spectral problem is already available. 
Negative coefficients /r®® are possible due to resonance effects. We mention that our 
approach could also be carried out with metallic inclusions {Eb with a negative real 
part), if one constructs resonators with a split ring structure as in [7j or |15j . 

The new feature in the present work is the network of thin wires. We have seen 
that this network contributes to the effective permittivity £®®. The formula fll.lll) 
for £®® is frequency independent, the relevant new contribution is explicitely given 
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as (and is not given through a cell problem). The wires do not create a neg¬ 

ative permittivity through some resonance effect, but merely through an averaging 
procedure: is the volume of the wires, eu,ri~‘^ is the permittivity in the wires. 

Nevertheless, let us emphasize that we observe here an effect that is more in¬ 
volved than some simple averaging: Only the connectedness of the wires across 
cells makes the effect possible (i.e.: the topology of the wires). Indeed, if T^ did not 
connect opposite sides, the test function could be constructed such that fl3.7p 
holds also in T)). In that case, the wire had no effect in the averaged law. 

A Bulk-resonance and the formula for /i®® 

In order to derive formula fll.lOp for /i®®, one has to calculate the F-averages of the 
solutions to the cell-problem of Lemma 12.21 We briefly sketch the arguments 
leading to fll.lOp . following [5]. The underlying concept of describing the cell- 
problem for with a bilinear form on a suitable Hilbert space has been used 
aleady in [7] (which was written earlier than [3]), but the useful concept of geometric 
averaging was only introduced with [3]. 

One considers the Hilbert space Xq := {u E H^iy) \ curlw = 0 in F \ Ey, § u = 
0} and the bilinear form bo{u,v) := /^{curlw • curlh -|- divM ■ divh}. The solutions 
of 02.191) are of the form = Cj -f Uj where Uj G Xq is determined by the 
variational equation (fc^ = cj^^oho) 

bQ{uj, v) — J Uj ■ V = k'^Eb J Cj ■ V \/v E Xq . 

The symmetric bilinear form 6o is coercive, it hence defines an operator Bq that 
has a compact self-adjoint resolvent on B^iY). The orthonormal eigenfunctions 
(7’n)neN to eigenvalues (A^)^^^ of Bq allow to express solutions as Uj = Cj^nf’n 
with Cj^n = Sbk’^i^n — Ebk‘^)~^ fy Cj ■ ifn- Definition 02.22p of the effective tensor /i®® 
provides fll.lOp . 
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